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these stars should have drastically reduced equatorial winds (mass loss) while 
on the main sequence, and as such should retain most of their angular momen- 
tum. Observations of both the Be/(B+Be) ratio and the blue-to- red supergiant 
ratio appear to have a metallicity dependence that may be caused by high rota- 
tional velocities. We have analyzed 39 archival Hubble Space Telescope Imaging 
Spectrograph (STIS), high resolution, ultraviolet spectra of O-type stars in the 
Magellanic Clouds to determine their projected rotational velocities Vsini. Our 
methodology is based on a previous study of the projected rotational velocities 
of Galactic O-type stars using International Ultraviolet Explorer (IUE) Short 
Wavelength Prime (SWP) Camera high dispersion spectra, which resulted in a 
catalog of Vsini values for 177 stars. Here we present complementary V smi 
values for 21 Large Magellanic Cloud and 22 Small Magellanic Cloud O-type 
stars based on STIS and IUE UV spectroscopy. The distribution of V sin i values 
for type stars in the Magellanic Clouds is compared to that of Galactic type 
stars. Despite the theoretical predictions and indirect observational evidence for 
high rotation, the type stars in the Magellanic Clouds do not appear to rotate 
faster than their Galactic counterparts. 



Subject headings: Stars: spectroscopic — stars: early-type — stars: rotation - 
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1. Introduction 

What effect does metallicity have on the mass loss and angular momentum loss of mas- 
sive stars? There is a large amount of secondary evidence, some of which is summarized 
below, suggesting that low metallicity should result in smaller angular momentum loss and 
subsequently higher rotational velocities for these type stars. However, there is an unfor- 
tunate lack of direct measurements of the equatorial rotational velocity V rot , or even of the 
projected rotational velocity Vsini for massive stars in a low Z environment. The Large 
and Small Magellanic Clouds (LMC & SMC) are nearby, contain many young massive ob- 
jects, and have low Z values of 0.007 and 0.004, respectively. As such they present an ideal 
laboratory to test the dependence of massive star rotation on metallicity. 

Maeder, Grebel, & Mermilliod (1999) found a significant increase in the number ratio of 
Be/ (B+Be) stars with decreasing metallicity among Galactic and MC clusters (see, however, 
the counterarguments presented by Keller et al. 1999). The Be phenomenon is closely linked 
to rapid rotation (Porter & Rivinius 2003). To explain this trend in Be star number ratio, 
the authors raise the possibility that stars forming in a low metallicity environment would 
have higher initial rotational velocities. Although no numerical models of the formation of 
massive stars at low Z have been produced, the authors tentatively present some possible 
origins of these higher initial rotational velocities. A lower Z environment would have lower 
dust content and fewer metallic ions present in star forming regions so that the magnetic field 
of the contracting central mass is less coupled to the surrounding region. Also the lifetime of 
the accretion disk is likely related to its opacity, which will decrease with lower metallicity. 
These factors would result in less angular momentum loss during star formation, leading to 
higher initial rotational velocities. The possibility of faster initial rotation at low metallicity 
resulting in more rotational mixing on the main sequence has also been suggested to resolve 
the nitrogen enrichment seen in SMC B- and A-supergiants (Venn et al. 1998). Maeder et 
al. (1999) state "Of course, direct observations of V'sinz in LMC and SMC clusters are very 
much needed in order to substantiate the above results." 

Observations of the number ratio of blue to red supergiants in nearby galaxies also 
show a sharp dependence upon metallicity (Humphreys & McElroy 1984). Stellar interior 
models which include rotation can account for this long-standing problem (Maeder & Meynet 
2001). The primary cause of this is the mild mixing just outside the core of a rotating star 
during the main sequence. This increases the amount of helium near and above the H- 
shell. The net effect of the larger He core and the mild main sequence mixing is to decrease 
the importance of the convective zone above the H-burning shell in the post-main sequence 
evolution. Because of the small polytropic index in the convective zone, the larger the size 
of the zone the smaller the radius of the star. Conversely, if its contribution is small, the 
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overall size of the star increases, leading to a redward position in the Hertzsprung-Russell 
(HR) diagram. 

Another key prediction of these new models is shown by comparing Figure 10 of Meynet 
& Maeder (2000) and Figure 3 of Maeder & Meynet (2001). These are plots of the evolu- 
tion of surface rotational velocities for massive stars in both Galactic and SMC metallicity, 
respectively. At low Z, the surface rotational velocities will remain almost constant during 
the main sequence evolution, contrary to that of massive stars at solar metallicity where 
the rotational velocities are predicted to rapidly decrease. Theoretically this results from 
both lower equatorial mass loss and the outward transport of angular momentum from the 
interior of the star by meridional circulation. Models of massive stars at both metallicities 
show a drastic decrease in rotational velocity after the terminal age main sequence. 

The idea that a low metallicity environment would create massive, rapidly rotating, 
helium stars has even been raised by theorists as a possible source of gamma-ray burst 
progenitors (MacFadyen & Woosley 1999). Here the outer layers of the star, and, in at 
least some cases, the mantle have too much angular momentum to fall freely inside the last 
stable orbit during core collapse. An accretion disk forms where the dissipation of rotational 
and gravitational energy will give rise to some sort of mass ejection and electromagnetic 
radiation, a "hypernova." This results in the prompt formation of a black hole through 
hyperaccretion. 

2. Projected Rotational Velocities from UV lines 

Stellar rotation is the dominant process in shaping the photospheric lines of O stars. 
The vast number of lines present in the ultraviolet spectra of O-type stars suggests that 
accurate projected rotational velocities could be found by studying these objects in the 
ultraviolet. These lines arise from high excitation transitions from deeper in the photosphere, 
and they are less subject to wind variability than optical lines. Also they are less likely to be 
contaminated by weak emission from circumstellar and nebular gas. The standard method of 
individual profile fitting in the UV is difficult because of line blending. A powerful alternative 
involves cross-correlation of a narrow-lined star with a test star, which will result in a cross- 
correlation function (ccf) that represents a "superfine" of the test star. This technique 
actually increases the S/N of the observed spectrum by in effect combining all the UV lines 
into one. The "superfine" will have an observed width related to the line width of the test 
star and that of the narrow-lined star (Figure 1). Fitting this "superfine" is the numerical 
equivalent of fitting each UV absorption feature and combining the resulting measurements. 
Our previously derived Vsmi values for 177 Galactic O stars were calculated by calibrating 
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the widths of ccfs derived from IUE SWP high dispersion spectra as a function of published 
V smi values from Conti & Ebbets (1977) in a common sample of objects (Penny 1996 = 
P96). Conti & Ebbets (1977) fit artificially broadened model line profiles (from the non- 
LTE model atmospheres of Auer & Mihalas 1972) of four visual absorption lines to determine 
their V sini values. The lines used were FEy A4340, He I A4388 and A4471, and He II A4541. 
Howarth et al. (1997) performed a cross-correlation study similar to ours with IUE high 
resolution spectra of 373 OB-type stars using the template star r Sco = HD 149438 (BO. 2 
IV). Although their calibration technique differed slightly from ours, their findings supported 
those from our study. 



3. Observations and Reductions 

Methods similar to those above can be utilized with observations of O-type stars made 
with the Space Telescope Imaging Spectrograph (STIS) E140M grating. The E140M grating 
produces in a single exposure an echelle spectrogram covering the wavelength interval 1150 — 
1700A, with an effective spectral resolving power of R = 46,000. This grating set up is 
complementary to the IUE SWP arrangement which was used so effectively to obtain an 
almost complete catalog of UV observations of Galactic O-type stars to V ~ 11. The 
Multimission Archive at Space Telescope (MAST 1 ) contains spectra of the 18 SMC and 6 
LMC O-type stars observed with this grating. There are also STIS E140M grating archival 
spectra of 12 Galactic O-type stars which we obtained in order to calibrate our methodology 
with the STIS observations. Although the E140M spectra are similar in many ways to 
the SWP high resolution spectra in our previous methodology, it is prudent to use these 
overlapping Galactic stars as an experimental checks of our calibration. There are also 
available in MAST IUE SWP high resolution spectra of 4 SMC and 16 LMC O-type stars; 
one of the latter is among those with archival STIS spectra. 

The NEWSIPS MXHI files of the above SWP spectra are spectrograms that cover 
the wavelength region of 1200 — 1900A. The archival STIS spectra are in extracted one 
dimensional FITS format, forming a near complete spectrogram covering the wavelength 
range 1150 - 1700A. We use a series of routines we have written in the Interactive Data 
Language 2 to further reduce the spectra (see details in Penny, Gies, & Bagnuolo 1997). One 
of the most important steps in our reduction is the placement of the data on a log A grid, in 
which each pixel step corresponds to a uniform velocity step of 1 pixel = 10 km s -1 . The 



1 URL: http://archive.stsci.edu/ 

2 IDL is a registered trademark of Research Systems, Inc. 
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Fig. I. — The cross-correlation functions and Gaussian fits for the SMC stars (a) NGC346- 
113 (OC6 V; V sini < 40 km s" 1 ), (b) AV 47 (08 III((f)); V sim = 76 km s" 1 ), (c) AV 423 
(09.5 V; Vsim = 186 km s" 1 ), (d) AV 80 (04-6n(f)p; l^sini = 324km s" 1 ). The bumpy 
nature of the ccf for AV 80 may be indicative of non-radial pulsations. 
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spectra were smoothed using a Gaussian transfer function (FWHM = 40 km s 1 ) and the 
interstellar features are removed. The spectra are then rectified to a unit continuum. 

Complete lists of spectra obtained, targets, and their spectral classifications are pre- 
sented in Tables 1, 2, & 3 for Galactic, SMC, and LMC targets, respectively. Spectral clas- 
sifications given with square brackets are from investigators other than Walborn (see Maiz- 
Apellaniz et al. 2004) and are taken from the compilations of Howarth & Prinja (1989), 
Conti & Ebbets (1977), Garmany, Conti, & Massey (1987), Conti, Garmany, & Massey 
(1986), Parker et al. (1992), Garmany, Massey, & Parker (1994), Massey et al. (1989), 
Massey, Parker, & Garmany (1989), and Massey et al. (1995). We make special note of the 
SMC star AV 423. This star is classified as 09.5 V by Massey et al. (1995). However, its UV 
spectrum, plotted in Figure 2, is inconsistent with that designation. Both Si IV AA1400 and 
C IV AA1550 doublets have P Cygni profiles, which makes this only the third star known in 
the SMC with Si IV emission. This feature in SMC is indicative of a supergiant classification. 
However, the N V doublet is extremely weak even in absorption. In addition to AV 423, we 
also plot AV 15 [06.5 11(f)], AV 83 [07 Iaf+], and NGC346-12 [O9.5-B0 V] for comparison 
in Figure 2. Whatever the revised type of AV 423, it is clearly not a dwarf star, and we will 
treat it as an evolved object below. 



4. Projected Rotational Velocities from Cross-Correlation Functions 

The spectra for each star are cross-correlated with a narrow-lined template star. We 
actually cross-correlate each test star with four different templates in order to get the best 
(usually highest) peak for fitting. These templates are HD 34078 (09.5 V; Vsini = 25 
km s" 1 ), HD 149438 (B0.2 IV; Vsini < 5 km s" 1 ), HD 54662 (06.5 V; Vsini = 85 km s" 1 ), 
and HD 57682 (09 IV; V sin i = 33 km s _1 ). We cross-correlate essentially the full wavelength 
range of the spectra (STIS: 1200 - 1700A; WE: 1200 - 1900A), setting to unity regions 
surrounding the strong P Cygni lines of N V A1240, Si IV A 1400, and C IV A 1550, since 
these features reflect wind and not photospheric rotational motion. The region around Lya 
is also set to unity. The spectra are padded with 1000 km s~ x of artificial continuum on 
both ends to include the entire observed range in the cross-correlation. The cross-correlation 
function (ccf) is the sum of the square of the differences between the test spectrum and the 
reference spectrum shifted in velocity from —1000 km s _1 to +1000 km s _1 at 10 km s _1 
intervals. The ccfs are rectified and inverted for convenience in fitting and plotting. Each 
ccf is then fitted with a Gaussian function, and V sin i values are derived from the Gaussian 
widths using the calibrations presented in P96 (HD 34078), Penny et al. (2001, HD 149438), 
Penny et al. (1997, HD 54662), and Penny (1997, HD 57682). We adopt 40 km s" 1 as our 
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Fig. 2. — Projected rotational velocities from STIS spectra versus those from IUE spectra 
for 12 Galactic O-type stars. The asterisk is plotted for the values for HD 13745, which may 
display time variations in line width (see §4). 
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minimum measurable Vsini, as this is the FWHM of the Gaussian transfer function used 
to smooth the data. 

We compared V sin % values from our previous study to those obtained using STIS ob- 
servations for those Galactic O-type stars in common. This is especially important since the 
wavelength range in the STIS spectra is smaller than that of the IUE observations, which 
might introduce systematic differences in our Vsini values. We plot in Figure 3 Vsini 
values from the STIS observations versus those from the IUE observations. We see excel- 
lent agreement, well within our estimated errors from the IUE data (±2 — 5%), with the 
exception of HD 13745. This star is noted in P96 as being a possible unresolved double-lined 
spectroscopic binary or it may be a short period pulsator (Koen & Eyer 2002), so a larger 
than expected difference in V sin i here is not unexpected. 

Our calculated projected rotational velocities are presented in Tables 2 and 3. Previous 
estimates of V sin i for the SMC stars in common by Walborn et al. (2000) are also shown in 
Table 2. We note that they describe their values as preliminary due to the blended nature of 
the He I + He II A4026 line from which their estimates are determined. Nevertheless, there 
is reasonable agreement between our results except in the cases of NGC346-324 and AV 47, 
which may be spectroscopic binaries observed at differing orbital phases. The SMC contains 
both the slowest and fastest rotators among our sample of Magellanic Cloud stars. Both 
NGC346-113 and AV 220 have narrow-lined spectra indicating a I/sinz below the minimum 
for which we can determine an accurate projected rotational velocity (< 40 km s _1 ). The 
fastest rotator is AV 80 that has a V sin i comparable to that found in some Be stars. Walborn 
et al. (2000) show that the He II A4686 feature is a double-peaked emission line in this star, 
which suggests that this star may have a disk-like circumstellar envelope, again like those 
found in Be stars. 



5. Results and Conclusions 

While our ultimate goal is to determine the true rotational velocities of O-type stars 
in differing environments, in general only the projected rotational velocities are measurable. 
However, we can make some statistical arguments. The inclination angles of the polar axes 
i of a group of stars should be randomly distributed. For each metallicity environment 
(Galactic, LMC, and SMC) we expect that the distribution of i should be the same, that is 
random. Thus, we can safely intercompare the rotational properties through an examination 
of their projected rotational velocities. We can compare the projected rotational velocity 
distributions of the different groups of star by calculating their cumulative distribution func- 
tion (cdf), as shown in Howarth et al. (1997) and Howarth & Smith (2001). The cdf gives 



- 9- 



Table 1. Galactic Targets 
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Table 2. Small Magellanic Cloud Targets 
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Fig. 3.— STIS spectra of AV 423 {offset by +6.0), AV 15 (06.5 11(f); offset by +4-0), AV 83 
(07 Iaf+; offset by +2.0), and AV 327 (09.5 II-Ibw). The current spectral classification of 
09.5 V for AV 423 is inconsistent with its UV spectrum. 
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Table 3. Large Magellanic Cloud Targets 
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the fraction of stars having a V sin i less than a specific upper value that ranges from zero 
to the maximum V sini observed. 

Before comparing the MC samples it is helpful to determine if the cdf of Galactic dwarf 
(IV- V) stars is different from that of Galactic evolved (I-III) stars. We plot in Figure 4 the 
CDF for both groups of stars based on the sample in P96. We see that the evolved sample 
has both fewer slow rotators and fewer fast rotators compared to the unevolved sample. The 
Kolmolgorov-Smirnov (K-S) test determines the probability that two samples are drawn from 
the same parent population. For the Galactic evolved and unevolved stars this probability is 
only 16%, which suggests that the differences may be significant. Howarth et al. (1997) came 
to the same conclusion in their study of projected rotational velocities of Galactic O-stars. 
They suggest that the evolved stars do indeed spin down as they grow in radius (removing the 
fastest rotators from the sample) and that there exists excess turbulent broadening among 
the supergiants (that removes the most narrow-lined stars from the sample). 

The differences in the distributions of projected rotational velocity between the un- 
evolved and evolved groups appears to support theoretical models that predict a spin-down 
during the core H- burning phases of evolution (Heger & Langer 2000; Meynet & Maeder 
2000). It is important to recognize that most stars with an O spectral type are probably 
found in the core H-burning phase, even among the more luminous stars (see the evolution- 
ary tracks of Schaller et al. 1992). Thus, we expect that the low Z, Magellanic Cloud O-stars 
are also in this phase where the rotational spin-down is predicted to be much less effective 
than we find in Galactic stars (Maeder & Meynet 2001). The low Z stars should retain an 
almost constant V rot during their core H-burning evolution. Therefore, we expected that the 
cumulative distribution function for the low Z stars would fall well below that of Galactic 
stars, because of the relative excess of rapid rotators. 

The existence of a systematic difference between the velocity distributions of Galactic 
unevolved and evolved O-stars (Fig. 4) indicates that we must compare like samples be- 
tween the Galactic and Magellanic Cloud stars in order to avoid spurious results because 
of differing proportions of unevolved and evolved stars in the groups. A large percentage 
of our Magellanic Cloud samples are objects with evolved luminosity classes (68% and 90% 
for the SMC and LMC samples, respectively). Note that we include the two SMC stars 
lacking a luminosity class, AV 80 and AV 220, among the evolved stars since Walborn et 
al. (2000) find them to have luminosities consistent with an advanced evolutionary state. 
We also include SMC star AV 423 among the evolved group because of the features in the 
UV spectrum indicating a high luminosity (see §3). We show in Figure 5 the cumulative 
distribution functions for the evolved O-stars in the Galaxy, SMC, and LMC. Surprisingly, 
the distributions for all three samples are very similar, with K-S probabilities of 44% and 




Fig. 4. — The cumulative distributions of projected rotational velocities for Galactic stars 
from Penny (1996). The dotted line represent those objects with luminosity classes IV and 
V, while the solid line corresponds to stars with luminosity classes I, II, and III. 
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42% that the LMC and SMC samples, respectively, are drawn from the same parent popu- 
lation as the Galactic sample. There is no visible trend of cumulative distribution functions 
with metallicity. The LMC sample has slightly faster V sin i values than that of the Galactic 
stars, but the SMC objects have slightly slower rotation rates than the other two samples. 
There are too few unevolved stars to make a reliable comparison, but the average projected 
rotational velocities in each group appear to be comparable: < V sini >= 131 ±93, 117±31, 
and 78 ± 45 km s _1 for the Galactic, LMC, and SMC unevolved stars, respectively. 

Taken at face value, our results do not support the prediction that the low Z stars of 
the Magellanic Clouds are rapid rotators. However, because of our relatively small sample 
sizes, we clearly need additional measurements to confirm the the result. Are low Z stars 
actually rotating at the same speeds as their Galactic counterparts? In order to determine 
conclusively the answer, a larger number of V sin i values for stars in these environments is 
needed. A previous study by Keller (2001) examined the rotation rates of 100 early B-type 
stars in the LMC. Compared with a similar sample of Galactic stars, the low metallicity stars 
show larger rotation rates. The largest initial mass of these early B-type stars is only 12 M Q 
and as such their MS evolution is not expected to be influenced by mass loss. However this is 
exactly the sort of large scale project which is necessary for the more massive O-type stars. 
A project with the FLAMES multi-object spectrograph on the VLT to obtain such spectra is 
currently ongoing. Additionally efforts are underway which focus on observing evolved (I-III 
classes) early type (03-05) stars in the LMC. As these stars should correspond to a later 
stage in the core H-burning evolution of the most massive stars, they represent a critical 
test of the current treatment of mass and angular momentum loss in low metallicity massive 
stars. 

Support for this project AR#09945 was provided in part by by NASA through a grant 
from the Space Telescope Science Institute, which is operated by the Association of Universi- 
ties for Research in Astronomy, Inc. under NASA contract NAS5-26555. The data presented 
in this paper were obtained from the Multimission Archive at the Space Telescope Science 
Institute (MAST). Support for MAST for non-HST data is provided by the NASA Office of 
Space Science via grant NAG5-7584 and by other grants and contracts. 
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